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Abstract 

Assuming the validity of grand canonical statistics, we study the properties of a spin- 
polarized Fermi gas in harmonic traps. Universal forms of Fermi temperature Tp, internal 
energy U and the specific heat per particle of the trapped Fermi gas are calculated as a 
function of particle number, and the results compared with those of infinite number particles. 
PACS numbers:03.75.Fi,05.30.Fk 



The ideal Fermi gas is an old and well-understood problem since the non-interaction 
Fermi gas is a good zeroth-order approximation for many familiar systems. Like the trapped 
degenerate atomic gases [1-3] that provide exciting opportunities for the manipulation and 
quantitative study of quantum statistical effects, the behavior of trapped Fermi gas also 
merits attention, both as a degenerate quantum system in its own right and as a possible 
precursor to a paired Fermi condensate at low temperature [4] though it perhaps is not as 
dramatic as the phase transition associated with bosons. 

The trapped atomic gases reported in Ref.[l-3] are dilute. The effects of predominantly 
short-range atom-atom interactions are therefore week. For dilute spin-polarized Fermi 
gases, the s-wave scattering amplitude vanishes due to the antisymmetry of the many- 
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Fermi wave function, and the p-wave scattering is small at low energy. Therefore the dilute 
Fermi gas is usually treated as an ideal Fermi gas. 

Unlike the usual system in most textbooks of statistical mechanics, e.g. [5], however, the 
number of trapped atoms in Ref.[l-3] is finite (about 10 9 ). For finite number of trapped 
fermions, the applicability of the usual thermodynamical calculations and the neglect of the 
ground energy are no longer available. For these relatively low numbers (compared to 10 23 ), 
the effects caused by the above approximation are nonvanishing. 

Harmonic traps provide a particularly simple realization of the confined Fermi system. 
In this trap, Butts and Rokhsar[6] calculated the spatial and momentum distributions 
for trapped fermions using Thomas-Fermi approximation. The spatial distribution of the 
trapped cloud provides an explicit visualization of a real-space "Fermi sea" and the mo- 
mentum distribution, unlike the spatial one, is isotropic. 

In this paper, we calculate the Fermi temperature, internal energy and the specific heat 
per particle of the finite number trapped Fermi gas. We indeed find marked differences 
from the usual treatments: a correction to the Fermi temperature, internal energy and the 
specific heat per particle of the finite number trapped fermions, a constraint to the number 
of trapped fermions for given trap frequencies. 

Consider N spin-polarized fermions of mass m moving in anisotropic oscillator poten- 
tials. The single-particle levels are familiar 

E N = E(n x ,n y ,n z ) = ^ fax^n* + ^), (1) 

i=x,y,z 

where uji(i = x,y,z)is the trap frequencies in the i direction, rij is non-negative integer. 
Note that each oscillator state is assumed to be filled with a single fermion, since only one 
spin orientation is confined by the magnetic trap. Under real experimental conditions of 
trapped atomic gases, the temperature is high on the scale of the trap level spacing, namely, 
UbT » h(Ji(i = x,y,z). Therefore, within the canonical ensemble, the partition function 

OO -i 

Q((3) = ^ e~P EN = e -f3(n x ui x +n y u> y +n z ui z ) _ J~J (2) 

N n x ,n y ,n z =0 i=x,y,z ^ e 
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of such a system without interactions can be expanded as follows: 

Q(P) ^ a 2 f3~ 3 + ai/T 2 + ao/T 1 + a_i + 0(/9) (3) 

with 

1 
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ai = ^02^+^ + ^) (4) 

a = -^a 2 (u^ + + + St^xtjj, + 3uj y uj z + 3uj z lu z ) 

a_i = - + —(u x /u y + u x /u z + Uy/u z +Uy/u x + u z /u x +u z /u;y), (5) 

where, /? = fe^T, &b is the Boltzmann constant, and the contribution of the ground state 
was singled away for special treatment. On the other hand, by making use of 

e-P E p{E)dE (6) 

iV -o 

the partition function can be calculated equivalently, where p(E) is the density of states, 
£q stands for the energy of the ground state. Comparing with eq.(3), it is proved that the 
density of states p(E) takes a form: 

p(E) = b + b 1 E + b 2 E 2 + .... (7) 

This expansion can't contain the terms E r with r being negative or non-integer, since they 
become zero after comparison with the direct calculation given by eq.(3). Substitution 
eq.(7) into eq.(6), one can easily find: 

Q{(3) = (b + he + b 2 el)(3' 1 + (h - 2b 2 eo)f3~ 2 + 26 2 /T 3 + ... (8) 

comparison of eq.(8) with the direct calculation(3) shows that 

bo = a - e ai - £q0 2 , h = aq + e a 2 , b 2 = 0.5a 2 (9) 
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There are two additional terms bo and 61 E in the density of states in comparison with the 
result of Butts [6], as we see in following, these additional terms will result in the shift of 
the thermodynamic quantities. 

Now, let us consider a system of N noninteracting fermions that the population N(Ei) 
of a state with energy Ei is given by the Fermi-Dirac distribution 

N ^ = &sha+i < 10 > 

Here, we set the statistical weights corresponding the state Ei, gi = 1. fi stands for the 
chemical potential, which is determined by the constraint that the total number of the 
particles in the system is N: 

00 

N = J2N(E t ). (11) 

i 

At zero temperature the Fermi-Dirac distribution factor is unity for energies less than the 
Fermi energy Ep = fi(T = 0,N), and zero otherwise. A straightforward integration of 
Eq.(ll) gives: 

E F + ^E F + 3^E F -^N' = (12) 

with N' = N + 1 62^0 + \pi £ 1 + bo£$. The terms with bo and b\ are corrections from finite 
number effects to Fermi energy. In the large number particle limit, Ep = {j^N') 1 ^ . This 
just is the result of Ref.[6,7]. Some words of caution are now in order. Fermi energy Ep 
is determined by eq.(12), i.e. Ep is a solution of eg. (12). From the mathematical point of 
view, the eg. (12) has a only real solution 

% 7 4 6 2 2 6 2 

in the situation of the number of trapped fermior satisfying 

N < N max , N max = [-(|) 3/2 + q ]b /3, (14) 
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where 

36? 60 b\ 36 &i 

P = -46| + 3 ^' 90 = 46|-^- (15) 

Otherwise, the eq.(12) have three real solution, but only the positive smallest one meets the 

requirements for Fermi energy from the physical point of view. We would like to point out 

that the results presented here are valid for k B T >> Twji(i = x,y,z). If this condition is 

broken, we should make use of a numerical method[8] to study the effects of finite particles. 

For finite temperature, neglecting the effect of zero-point energy, the eq.(ll) gives: 

N = b 2 k s B T 3 f 3 (z) + hklT 2 /^) + b k B Th(z) (16) 

where f n (z) stands for the Fermi integral, which is given that 

1 poo T n—1 poo 

ff z )= —Z dx,T(n)= e~ x x n - x dx. 

JnK ' r(n) Jo z~ 1 x n + 1 v ' Jo 

z is a fugacity. The last two terms in eq.(16) are due to the effects of finite particles. If one 

introduces a temperature 

t° f = ^(r^) 173 - (17) 

which denotes the Fermi temperature of infinite number fermions trapped in anisotropic 
oscillator potentials [6], i.e. k B T F = fi. Then the Fermi temperature Tp defined as the same 
as in most textbooks for finite number fermions in the case of £0 = takes a form 

T -r°ri lbl 1 2b ° 1 jjjggO i (-1^ 

lF lF[L 3b 2 k B T F h(z F ) 3b 2 (klT F ff 3 (z F ) U 

where The last two terms in eq.(18) are correction to Fermi energy from 

the finite number effect, it shows that the Fermi temperature for finite trapped fermions is 
lower than that in the case [7] of infinite trapped fermions. 

Generally speaking, the specific heat is more interest from the experimental point of 
view, since the low-temperature behavior of specific heat c is generally treated as the hall- 
mark of onset of phase transition . It is well known that the specific heat can be derived 
from the internal energy, which is expressed by 

pu=pr jgcgL dE (19) 

7 exp[/3(E - e )\ + 1 V ' 
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Here we dropped the ground state energy, which is less than 2eq and remains constant. 
Substituting eq.(7) into eq.(19), one easily finds 

U= b fj,{z)+ b ^h{z)+ b p 2 {z), (20) 

the first term remains unchanged in the limit of very large number particles, while the last 
two terms close to zero. The specific heat per particle of the trapped fermions is given 
that [6] 

r 1 9U, 

c=[^k, (21) 

straightforward calculation gives 

c = ^{4b 2 k 4 B T 3 U(z) + ^kl^hiz) + 2b k 2 B Tf 2 (z) - 

3 ,4^ /3/2(^/3^) 3 3 2 /3/2(^/2^) 3 2 f 3 /2(z)fi(z) 

2 hkBT ~Mi) 2 hkBT 2 b ° kBT f 1/2 (z) } (22) 

It is evident that the specific heat per particle has increased due to the finite particle effects. 
The results are illustrated in Fig.l. As state in most textbooks of statistical mechanics, for 
ideal Fermi gas in a 3D box, the specific heat increases proportionally with temperature 
T for T << Tp, whereas it closes to a constant 1.5/cb for T — > 00. For trapped fermions, 
however, the specific heat increases predominantly with T for T << Tp, and close to 3&b 
with T — > 00. Physically, the difference between the case of 3D box and that of trapped 
fermions is the confinment. This leads to the difference in their energy levels. Consequently, 
different energy levels would result in different dependence of specific heat on temperature. 

In conclusion, we have discussed the behaviors of fermions trapped in an anisotropic 
oscillator potentials, it was shown that corrections due to the effect of finite fermions and 
the ground state energy are small, whether they are observable depends on the experimental 
parameters. We believe the results will be useful to them who is working experimentally in 
this area. 
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Fig.l The specific heat per particle vs. temperature. The dots show the situation 
for a very large number of particles(N = 10 23 ), and the dashed line that for a smaller 
number(10 8 ). The trap frequencies are to x = 500Hz,u> y = 600Hz,lj z = 800Hz. 
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